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I. INTRODUCTION 

Superposition and entanglement are properties that 
distinguish quantum mechanics from any classical the¬ 
ory. In an entangled quantum system there are states 
that cannot be factorized into product of states of the 
subsystems. Entanglement has became of interest since 
it is a key property in quantum information, cryptogra¬ 
phy and quantum computation [lj, Q. There exist several 
sources of entangled quantum systems. For instance, in 
solid state physics, quantum optics and also atoms in 
cavity quantum electrodynamics. In the literature many 
ways were proposed to generate entangled states in sys¬ 
tems of two-level atoms interacting with a bosonic field. 
See for example SH3- 

Researches on atomic radiative processes have proved 
to be of major importance in quantum optics. Following 
the seminal work of Dicke Q, there is an extensive liter¬ 
ature on resonant interaction of atoms [sHU- Recent 
examples, of great current interest, concern investiga¬ 
tions of radiative processes involving entanglement. Let 
us briefly discuss some of such important works. For in¬ 
stance, in Ref. [3 the authors investigate the properties 
of emission from two entangled atoms coupled with an 
electromagnetic field in unbounded space. In Ref. [3 the 
authors study the radiative processes of entangled two- 
level atoms coupled individually to two spatially sepa¬ 
rated cavities. Ref. [3 presents an interesting scheme to 
realize a highly controlled and selective atom-field inter¬ 
action in cavity quantum electrodynamics system. Non¬ 
linear optical processes in chains of ions in an entangled 
state were investigated in Ref. Q3- See also Ref. [3 • We 
mention that radiative processes of maximally entangled 
states are quite different from the non-entangled states. 
For example, in an entangled two-atom system there are 
superradiant and subradiant states. These states have 
very different spontaneous decay rates. The symmetric 
state decays with an enhanced whereas the antisymmet- 
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ric states decays with a reduced spontaneous emission 
rate fliti ]. 

By using first-order perturbation theory it can be 
shown that the transition rate of an atom interacting 
with an electromagnetic field in the vacuum state is 
given by the Fourier transform of the positive frequency 
Wightman function evaluated on the world line of the 
atom S!3- The heuristic picture is then that the atom 
is forced to radiate by the vacuum field fluctuations. On 
the other hand, following the discussion by Ackerhalt, 
Knight and Eberly [3 it is possible to interpret sponta¬ 
neous decay as a radiation-reaction effect. As discussed 
by Milonni, both effects, vacuum fluctuations and radi¬ 
ation reaction, depend on a particular ordering chosen 
for commuting atomic and field operators [HI] ■ Hence, it 
turns out that the contributions of vacuum fluctuations 
and radiation reaction can to a large extent be chosen 
arbitrarily. On the other hand, Dalibard, Dupont-Roc, 
and Cohen-Tannoudji argued that there exists a preferred 
operator ordering: when one chooses a symmetric order¬ 
ing the distinct contributions of vacuum fluctuations and 
radiation reaction to the rate of change of an atomic ob¬ 
servable are separately Hermitian and hence they possess 
an independent physical meaning [HI [HI . This became 
known as the DDC formalism. Within such a formalism 
the interplay between vacuum fluctuations and radiation 
reaction can be considered to maintain the stability of the 
atom in its ground state. We remark that this formalism 
was employed in many situations [Hi - [H| - The essential 
picture that emerges from such analysis is the following. 
The fluctuations of the quantum field act on the atoms, 
causing their polarizations; the atoms then fluctuate and 
disturb the field, which in turn reacts back on the atoms. 
It is clear that the proper comprehension of the interplay 
between both effects is crucial in the studies of resonant 
interaction of atoms. 

The dynamics of a small system, entangled or not, cou¬ 
pled to a reservoir can be discussed using different for¬ 
malisms mi. The key point is that the reservoir is gen¬ 
erally described by an infinitely many degrees-of-freedom 
formalism. The standard formalism for the evaluations 
of time evolution and correlation properties of collective 
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atomic systems is the traditional master equation ap¬ 
proach. Within this approach it is possible to quantify 
the degree of entanglement of a particular state. Here 
since we are interested in understanding the mechanism 
responsible for supporting entanglement in radiative pro¬ 
cesses involving atoms, we choose a different route and 
use the DDC formalism aforementioned. For an interest¬ 
ing investigation regarding boundary effects on quantum 
entanglement, see for instance [32| . 

We choose to go focusing our attention in applying this 
method to quantum entanglement of atoms. Being more 
specific, we are interested to study the contributions of 
vacuum fluctuations and radiation reaction to the stabil¬ 
ity or decay of an entangled state. We remark that in 
the present paper we do not consider the rotating-wave 
approximation. The organization of the paper is as fol¬ 
lows. In Section II we discuss the Hamiltonian describing 
a system of bound atoms and radiation field. In Sec¬ 
tion III we calculate the rate of variation of the atomic 
energy in vacuum when both atoms are at rest. Con¬ 
clusions and final remarks are given in section IV. We 
briefly review the discussion on relevant electromagnetic 
correlation functions in the Appendix. In this paper we 
use units H = c = ks = 1- We are using the Minkowski 
metric r] a p = — l,a = (3 = 1,2,3, rj a p = l,a = /3 = 0 
and r] a 0 = 0, a ^ (3. 


II. TWO IDENTICAL ATOMS COUPLED WITH 
AN ELECTROMAGNETIC FIELD 


Let us consider two globally neutral systems of charges 
localized at fixed distances iq and r 2 interacting with a 
quantum electromagnetic field. In what follows we will 
interpret such systems as identical atoms. We are work¬ 
ing in a four-dimensional Minkowski space-time. The 
Hamiltonian in the Coulomb gauge describing this sys¬ 
tem in the dipole approximation reads fifl 
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where Vc* ul (Vc„ u1 ) is the respectively Coulomb energy 
of the system of charges and dip is the electrostatic 
interaction energy between the atomic electric dipole mo¬ 
ments pi and p 2 of the systems of charges - the dipole- 
dipole interaction energy. The sums are over all charges 
comprising each of the systems. The last part is the 
Hamiltonian of the free electromagnetic field, in which 
a kA> a k.A represent the usual creation and annihilation 
operators of the electromagnetic field. In the above we 


have neglected the rest mass energy of the systems of 
charges and the zero-point energy associated with the 
quantum field (the energy of the longitudinal field can be 
considered as a correction to the mechanical rest mass). 
One expects the dipole approximation to be accurate 
whenever the coupling between particles and radiation 
mostly encompasses the modes whose wavelength is much 
larger than the typical size of the system. The above ex¬ 
pression is known as the minimal-coupling form of the 
Hamiltonian. 

The above Hamiltonian can be converted to a more 
convenient and plain form under the Power-Zienau- 
Woolley transformation [3 [35|. One gets 
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( 2 ) 


The dipole-dipole interaction energy has been compen¬ 
sated for by dipole terms which also give rise to dipole 
self-energies. The latter have been neglected in the above 
expression since it does not contribute to the coupling be¬ 
tween the system of charges and the quantum held. In 
fact, one can show that the contributions of cross terms 
related with the transverse polarizations of the systems of 
charges cancel the contribution coming from L^ 2 dip 33j. 
The total electric held in the expression above is to be 
interpreted as an electric displacement. One should un¬ 
derstand that the coupling term p^ • E(iqi) contains the 
interaction of the system of charges A with the trans¬ 
verse electric (displacement) held as well as the longi¬ 
tudinal held generated by the other system of charges 
at r a- The consequence of such considerations is that 
there are no more instantaneous electrostatic interac¬ 
tion terms between the systems of charges; all interac¬ 
tions now are realized through the quantum electromag¬ 
netic helds which propagate with velocity 1. This form 
is suitable for describing retarded dipole-dipole interac¬ 
tions between the atoms. Expression © is known as 
the multipolar-coupling form of the Hamiltonian. This is 
the representation that we shall consider throughout this 
paper. 

Hereafter we will be working within the Heisenberg pic¬ 
ture. As mentioned above, we assume that the globally 
neutral systems of charges describe two identical two- 
level atoms at rest located at iq and r 2 , respectively. We 
denote by |g,) and |e») the ground and excited state of 
the z-th atom, with energies =fwo/ 2, respectively. We 
employ this notation in order to allow the possibility of 
a straighforward generalization of our results to the case 
of non-identical atoms. The time evolution of the total 
system is to be taken with respect to the proper time r of 
the atoms. The important fact to be noticed is that the 
purely atomic parts that appear in equation © describe 
the free Hamiltonian of the two atoms labelled as Ha(t). 
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In the Dicke notation [§J, one has 

Ha(t) = ^ [{Sf ( r ) ® i ) + (i ® Si(r))] , ( 3 ) 

where S* = \e a ){e a \ - \g a ){g a \ is the energy operator 
of the a-th atom. The space of the two-atom system 
is spanned by four product stationary states which are 
eigenstates of Ha with respective energies 

E gg = -uj 0 \gg) = \gi)\g 2 ), 

E ge = 0 \ge) = \gi)\e 2 ), 

E eg = 0 \eg) = |ei)|g 2 ), 

E ee = w 0 |ee) = |ei)|e 2 ), (4) 

where a tensor product is implicit. Instead of working 
with this product-state basis, we can conveniently choose 
the Bell state basis. In terms of the product states, one 
has: 

\n ± ) = ^=(\g 1 )\e 2 )±\e 1 )\g 2 )) 

\$ ± ) = -±=(\g 1 )\g 2 )±\e 1 )\e 2 )). (5) 

The Bell states are known as the four maximally entan¬ 
gled two-qubit Bell states, and they form a convenient 
basis of the two-qubit space. Nevertheless, we remark 
that the Bell states ji^) are not eigenstates of the atomic 
Hamiltonian Ha', in turn, in view of the degeneracy as¬ 
sociated with the eigenstates \ge) and | eg), any linear 
combination of these degenerate eigenstates is also an 
eigenstate of the atomic Hamiltonian corresponding to 
the same energy eigenvalue. Therefore, the Bell states 
l^) are eigenstates of Ha- 
The stationary trajectories guarantees that the undis¬ 
turbed atomic system has stationary states. In the fol¬ 
lowing we will be interested in the evolution of the quan¬ 
tum electromagnetic field with respect to the proper time 
t. In this way, the Heisenberg equations of motions tells 
us that 

= E / (4^3 Wka k,A(*( r )) a k,A(t(T)) (6) 

Finally, from equation © one can easily identify the 
Hamiltonian which describes the interaction between the 
atoms and the field in the dipole approximation: 

Hi(t) = —Mi(r) • E(a:i(r)) - (r) ■ E(x 2 (r)). (7) 

The dipole moment operator is endowed of only off- 
diagonal elements and hence can be written as 

Mi(r) = {gi\f J 'i\ e i)\gi)( e i\ + H.c. = /x [Sf (t) + S~ (r)] , 

( 8 ) 

(no summation over repeated indices) where we have as¬ 
sumed a proper choice of phases which allows the dipole 
matrix elements to be real and we have defined the 
dipole raising and lowering operators as Sf = \(-i)(gi\ 


and S~ = |(?i)(ej|, respectively. In addition, since the 
atoms are presumed to be identical and similarly ori¬ 
ented, (gi\fj,i\ei) is independent of the index i and de¬ 
noted simply by fi. Incidentally, suppose that our atoms 
are spinless one-electron systems. Hence MaM = er a (r), 
where e is the electron charge and r a (r) is the position 
operator of the atom a. 

As discussed above, there is an extensive literature in¬ 
vestigating the role of vacuum fluctuations and radiation 
reaction in radiative processes of atoms. Here our main 
intention is to identify the distinct contributions of quan¬ 
tum field vacuum fluctuations and radiation reaction to 
the entanglement dynamics of the atoms. For such pur¬ 
poses, let us briefly discuss the DDC formalism. The first 
step is to set up the the Heisenberg equations of motion 
for the dynamical variables of the atom and the field with 
respect to r. Such expressions can be derived from the 
total Hamiltonian H(t) = Ha(t) + Hf{t) + Hj{t). Af¬ 
terwards one splits the solutions of the equations of mo¬ 
tion in two parts, namely: The free part, which is present 
even in the absence of the coupling; and the source part, 
which is caused by the interaction between atoms and 
field. Therefore one can write, for the atomic observables 
and the dynamical variables of the field, respectively: 

S z a (r)=S^(r) + S z a ’ s (r), 

Ok ,a(£(t)) = a[ x (t(T)) + < iA (*(t)). (9) 

a = 1,2. Since one can construct from the annihilation 
and creation field operators the free and source part of 
the quantum electric field, one also has 

E(f(r), x(t)) = E^(t(r),x(r)) + E s (t(r), x(r)). (10) 

Nevertheless, such a procedure introduces an ambiguity 
of operator ordering. The source part of the field acquires 
contributions of atomic observables during its time evolu¬ 
tion which implies that the fundamental aspect that all 
atomic observables commute with field variables is not 
preserved in time for E-^ and E s separately. Because the 
above expression contains products of atomic and field 
operators, one must choose an operator ordering when 
discussing the effects of E-^ and E s separately. Differ¬ 
ent operator orderings will yield different interpretations 
regarding the roles played by vacuum fluctuations and 
radiation reaction. As pointed out in the Introduction, 
only the adoption of a symmetric ordering of atomic and 
field variables permits both contributions to be Herrni- 
tian. In this way, the effects of such phenomena can 
posses an independent physical meaning. Therefore, fol¬ 
lowing this prescription, one can clearly separate the con¬ 
tribution of the vacuum fluctuations from the radiation- 
reaction effects in the evolution of the atoms’ energies, 
which are given by the expectation value of Ha, which 
in turn is given by equation J3]). Furthermore, in a per¬ 
turbative treatment, we take into account only terms up 
to order e 2 . We also consider an averaging over the 
field degrees of freedom by taking vacuum expectation 
values. In turn, since we are interested in the evolu¬ 
tion of expectation values of atomic observables, we also 
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take the expectation value of the above expressions in 
a state |w), which can be one of the states given by 
equation 0 or one can also consider the Bell states (0. 
The details of the calculations can be found elsewhere, 
see for instance {26]. Therefore, employing the notation 
(('“)) = (0, w|(- • • )|0, w), one has the pivotal results con¬ 
cerning the vacuum-fluctuation contribution 

r - 2 

I dr' ^2 Dij(x a (T),x b (T')) 

T ° o,6=l 

x ^A«(t,t'), (11) 

where 

A a&( T > r ') = a,& = l,2, (12) 

is the linear susceptibility of the two-atom system in the 
state |w) and 



Dij(x a (T),x b (T')) = (0| {e{ (,x a (r)),E / (x&(t'))}|0), 

(13) 

is the Hadamard’s elementary function. Since we are 
dealing with free fields, one can employ the results de¬ 
rived in the Appendix [A] On the other hand, for the 
radiation-reaction contribution, one has: 

/- J-=\\ dT ' A ij(x a {r),x b {T')) 

\ / RR Jto a,6=1 

X (ry), (14) 

where 

D ab( T ’ T ') = ( w l{/ i a / (' r )iM6 / ( r, )}|w), a,6= 1,2, (15) 

is the symmetric correlation function of the two-atom 
system in the state |w) and 

A l;j {x a {T),x b {T')) = <0| [e{( x a {T)),E f 3 (x b (r'))} |0), (16) 

is the Pauli-Jordan function which, for the same reason as 
above, can be evaluated from AppendixjA] Note that A®^ 
and D'J h characterize only the two-atom system itself, 
unlike the statistical functions of the field which have to 
be evaluated along the trajectory of the atoms. We see 
from equations CD and (fl4l) that the rate of variation of 
the energy of the two-atom system presents contributions 
from the isolated atoms and also contributions due to 
cross correlations between the atoms mediated by the 
field. This interference is a consequence of the interaction 
of each atom with the field. 

Using that ^ (r) = e lffA ' r /r( t ’l(0)e _ * ffAT , and 

the completeness relation for the two-atom states 


l a;, )( a;, | = 1; one can show that the statistical func¬ 
tions for the two-atom system can be put in the following 
explicit forms 
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where we have defined 


■m 1 >v) = (is) 


In this paper we are primarily interested in transitions 
from entangled states [currently represented by the Bell 
states 0 ] to separable states [given by equation 0 ] , or 
the inverse. For the purpose of investigating the fine de¬ 
tails regarding the degradation of entanglement between 
the atoms as a spontaneous emission phenomenon, as¬ 
sume that the atoms were initially prepared in one of the 
Bell states |r2 ± ). With this regard, equations 0 and 0 
state that the only allowed transitions are [lu) —> \gg), 
with Aw = w — u/ = wo > 0 and lU 1 * 1 ) —> |ee), with 
Aw = uj — oj' = —wo < 0. On the other hand, suppose 
that one is interested to address the generation of entan¬ 
glement. Assume that the atoms were initially prepared 
in the excited state |ee). For the Bell states |$ ± ) one 
gets Aw = 0, which implies that ( dHA/dt) = 0. Hence 
it is not possible to generate such Bell states out of the 
excited atoms. One then must consider the decay rate to 
one of the Bell states |f2 =l= ), obtaining Aw = wq > 0. 


III. RATE OF VARIATION OF THE ATOMIC 
ENERGY IN VACUUM FOR STATIC ATOMS 

In this Section we will consider in detail the rate of vari¬ 
ation of atomic energy for an inertial two-atom system us¬ 
ing the technique developed above. We assume that the 
atoms are separated by the space-time interval Ax M = 
(x(n) - x(t 2 )Y = (n - r 2 ,x 0l - x 02 ) = (ti - r 2 , Ax 0 ). 
The contributions m and m to spontaneous emis¬ 
sion can be evaluated from the results discussed in the 
Appendix [A] Being more specific, since x = x(t) and 
x' = x(t') we insert in such expressions the statistical 
functions of the two atom system, given by equations ED 
and ED, and the electromagnetic-field statistical func¬ 
tions given by (IA3I) and (1A6I) . Initially let us present 
the contributions coming from the vacuum fluctuations. 
These are given below, with u = r — t'\ 
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At 
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2 (Ax)i(Ax)j - Sij [|Ax| 2 + (u - ie) 2 ] 

2(Ax)j(Ax)j - [|Ax| 2 + (u + ie) 2 ] 


it .r) 

i / VF u/ a,6=1 -'-At 


[(u — ie) 2 — | Ax| 2 ] -= 


[(« + ie) 2 — | Ax| 2 ]" 


( 20 ) 


where At = t — tq, Ax = Xo Q — Xo b . By invoking the method of residues one may compute all the relevant integrals. 
For e —> 0 one gets 


' dJ t) 


47T 

VF u/ a,b— 1 

x ^0(Aw) — 6 {— Aw)j 

( ( Ax )i( Ax )t 


l Ax l 

( 3(Ax)j(Ax)j 


l l A > 


2 % 


1 — Aw|Ax| cot(Aw|Ax|) 


l Ax l : 


V l A > 


- l A 4 


( 21 ) 


where we have taken At —> oo. Let us note from the above contributions that vacuum fluctuations stimulate excitation 
(( dHA/dt) > 0) as well as deexcitation (( cLHA/dt) < 0) of the two-atom system. This is reminiscent from the fact 
that for an atom interacting with quantized radiation, one has that stimulated excitation and deexcitation have equal 
Einstein B coeffcients, which implies that vacuum fluctuations tend to excite an atom in the ground state as well as 
deexcitate an atom in the excited state. In our context, following the discussion at the end of Section ITU one realizes 
that vacuum fluctuations can promote the degradation of entanglement of atoms as well as activate the generation of 
entanglement between the atoms. 

Now let us evaluate the radiation-reaction contributions. One has, with u = t — t'\ 


l~^~) = M ab( u ^') Au} [ dMsin(Awu) 

\ / 1 


RR 


uj r a,6=1 
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| Ax| 2 
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By employing standard Dirac delta-function identities we get: 

/cLHaV 1 ^ ^ Aw sin(Aw|Ax|) 

\ dr J Air A -' 


47T 

RR w 0,6=1 


| Ax | ( 0 ( Aw ) + 0(-Aw)).M l o J b (w,w') 
3< ~ A |Axp X ^ “ i 1 ~ Aw l Ax l cot(Aw|Ax|)) ■ 


| Ax 

(Ax)i(Ax)j 
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l Ax l 2 


- 5^ (Aw) 2 


( 22 ) 


(23) 


Note that our results explicitly show contributions from 
a classical phenomenon (radiation-reaction effect) to the 
entanglement dynamics of two identical atoms. Adding 
up the above terms leads to the main conclusions of this 
work. First of all, we note the appearance of spatial oscil¬ 
lations in the rate of variation of atomic energy. In turn, 
we observe that for large times the difference between 


vacuum fluctuations and radiation reaction arises when 
considering excitation processes. In the perspective of 
spontaneous emission, both give identical contributions. 
Thus only emission processes are allowed for asymptotic 
times. Being more specific, the Bell states l^) are stable 
for absorption processes but through spontaneous emis¬ 
sion the atoms can disentangle and decay to the ground 
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state. On the other hand, the transition \gg) —> |fi ± ) 
is suppressed for At —> oo. This can be seen as a di¬ 
rect consequence of the known result that an atom in its 
ground state is stable due to the balance between vacuum 
fluctuations and radiation reaction. In order to perturb 
such a balance for At —> oo one must consider general 
non-inertial trajectories, such as uniformly accelerated 
motion [26] . Notwithstanding one can notice that entan¬ 
glement between the atoms can still be created, but only 
via emission processes; for instance, one must initially 
prepare the atoms in the excited state |ee) in order to 
detect the formation of entangled states for At —> oo. 
In summary, the importance of such results lies in the 
fact that by considering the cooperation between vacuum 
fluctuations and radiation-reaction effect, one could give 
a detailed mechanism that accounts for the generation as 
well as the degradation of entanglement in the resonant 
interaction of atoms. 

For completeness let us present the total rate of change 
of the atomic energy for At —> oo assuming that 
the atoms were initially prepared in one of the Bell 
states If^). This is obtained by adding the vacuum- 
fluctuations and radiation-reaction terms. For simplicity 
consider that the atoms are polarized along the X 3 di¬ 
rection. We also assume the special case in which the 
atoms are separated by the space-time interval Ax^ = 
(x(ti) - x(t 2 )Y = (n -t 2 ,0,0, (xi ) 3 - (x 2 ) 3 ). The total 
rate of change of the atomic energy is zero for Aw < 0. 
As discussed above the distinct contributions for sponta¬ 
neous excitation cancel each other in this situation, leav¬ 
ing only the terms for which Aw > 0. In addition, from 
equation ( J 4 |) and the aforementioned assumption on the 
dipole matrix elements one has that 

M\{(n±,g) = ^f 
Mg(n ± ,g) = ^f 

M%(fl ± ,g)=M 21 (il ± ,g)=±^. (24) 

Hence the explicit result for the rate of variation of energy 
of the two-atom system is given by: 


'dH A \ 
cLt / 

1 / t 


-^[l±/(w 0 Az)](p, 3 ) 2 , (25) 


where the (minus) plus sign refers to the (anti)symmetric 
Bell state (|H - )) |H + ) and Az = (x^ — (x 2 ) 3 - Also, we 
have defined the function: 


. 3 sin x / \ . . 

f(x) = -g—^1 — xcotxj. (26) 


Such a function quantifies the influence of the cross cor¬ 
relations on the entanglement between the atoms. Some 
special values are given by 


/(titt) 


3 

(mr ) 2 


(n is a positive integer) and 


/((n+l/2 W = (n+1 3 / 2 ) , j3 (-ir. (28) 

The behavior of this function is depicted in figure [T| 
Note the great oscillatory regime for short distances 
between the atoms. In fact, the plot shows that for 
|Az| -C 1/wo cross correlations are more important for 
the rate dH A /dx in comparison with the case in which the 
distance between the atoms is very large (|Az| 1/wo). 
It means that the cross correlations generate a construc¬ 
tive interference when the atoms are very near each other 
and these interference terms vanish for large spatial sepa¬ 
rations between entangled atoms. This is a similar result 
as described in references [TlL HB|. On the other hand, 
equation ( 125 T ) shows that, for small distances between the 
atoms, the symmetric Bell state |H + ) decays with an en¬ 
hanced emission (superradiant) rate as compared with 
terms of isolated atoms, whereas the antisymmetric Bell 
state |fl _ ) presents the opposite behavior (subradiant). 
This clearly indicates that the antisymmetric entangled 
state decays in a softened rate. In addition, note that 
for |Az| <C 1/wo the antisymmetric Bell state presents 
a nearly complete inhibition of the decay to the ground 
state due to destructive interference of cross correlations 
between the atoms. Such results confirm the assertions 
discussed in the Introduction. 



FIG. 1. The quantity f(x) as a function of x = wqAz. 


(-l ) n+1 


(27) 
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IV. CONCLUSIONS AND PERSPECTIVES 

Within first-order perturbation theory the transition 
rate of an atom interacting with an electromagnetic field 
in the vacuum state is given by the Fourier transform of 
the positive frequency Wightman function evaluated on 
the world line of the atom. Within this picture the atom 
is forced to radiate by the vacuum field fluctuations. On 
the other hand, it is possible to interpret spontaneous de¬ 
cay as a radiation-reaction effect. Both effects, vacuum 
fluctuations and radiation reaction, depend on a partic¬ 
ular ordering chosen for commuting atomic and field op¬ 
erators. On the other hand, as discussed above there 
exists a preferred operator ordering: when one chooses 
a symmetric ordering the distinct contributions of vac¬ 
uum fluctuations and radiation reaction to the rate of 
change of an atomic observable are separately Hermitian 
and hence they possess an independent physical mean¬ 
ing [24], [25| . Within such a formalism the interplay be¬ 
tween vacuum fluctuations and radiation reaction can be 
considered to maintain the stability of the atom in its 
ground state. Here in this paper we studied the contri¬ 
butions of vacuum fluctuations and radiation reaction in 
radiative processes of entangled atoms. We found that, 
at asymptotic times, with respect to excitation processes 
such contributions cancel each other whereas both mech¬ 
anisms are responsible for the decay of the entangled 
state. One possible way to unbalance the contributions 
from vacuum fluctuations and radiation reaction in ab¬ 
sorption processes is to consider situations of finite ob¬ 
servation time intervals. The main conclusion one can 
draw from such results is that when considering radia¬ 
tive processes of atoms the machinery underlying entan¬ 
glement can be envisaged as an interplay between clas¬ 
sical concepts (radiation-reaction effect) and quantum- 
mechanical phenomena (vacuum fluctuations). We be¬ 
lieve that the results presented in this paper may have 
an impact in the studies of radiative process of atoms. 
For instance, recently the subject of van der Waals in¬ 
teraction between atoms was taken up in reference (37] 
in which the authors demonstrate that such an inter¬ 
action presents both temporal and spatial oscillations. 
In the experimental side, the angular dependence of the 
resonant dipole-dipole interaction between two individ¬ 
ual Rydberg atoms with controlled relative positions was 
measured [381 ]. A framework in which vacuum fluctua¬ 
tions and radiation-reaction effect have been clearly iden¬ 
tify and quantitatively analyzed may contribute to an 
accurate and deeper understanding of such results. 

Conducting plates or cavities can modify the rate of 
spontaneous emission of excited atoms |39[. Since the 
presence of boundaries affect the vacuum fluctuations 
of the quantum held one should expect that the tran¬ 
sition rates of atoms are modified in this situation fio} . 
Hence, for entangled atoms, it is interesting to ask how 


the rate of variation of the atomic energy is modified 
in a situation where translational invariance is broken. 
Recently, boundary effects on entanglement were investi¬ 
gated assuming a framework where the spontaneous de¬ 
cay of atoms is only attributed to vacuum-fluctuation 
effects f3(|. A natural extension of our results is to dis¬ 
cuss the mean life of an entangled state in the presence of 
boundaries within the formalism employed in this work. 
This is under investigation by the authors. 
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Appendix A: Correlation functions of the 
electromagnetic fields in the vacuum 


In this appendix we shall concisely discuss free elec¬ 
tromagnetic correlation functions. We consider a four 
dimensional unbounded Minkowski space-time. In ad¬ 
dition, as previously mentioned we will work in the 
Coulomb gauge. Here we will not discuss the quantiza¬ 
tion of the electromagnetic field in detail. For a thorough 
analysis of quantum electrodynamics in the Coulomb 
gauge, we refer the reader the reference 0. In turn, 
details concerning the correlations of the electric field for 
two static space-time points can be found in the com¬ 
prehensive account by Takagi [41} and related references 
therein. Instead, we limit ourselves to present all relevant 
electromagnetic correlation functions in the vacuum em¬ 
ployed in the text. 

Let us start with the correlation functions of the elec¬ 
tric field. Through standard procedures, one may eval¬ 
uate the Pauli-Jordan function of the electric field. It is 
given by 


A«(i,i') = {0|[£‘W,£Vi')]|0) 

A (Ai) 


where 

D(t,x) = 
1 


47t|x| . 
lim 


<K|x| ~t) -£(|x| +t) 


(27t) 2 |x| e-s-o 


+ (|x|-i) 2 e 2 + (|x|+t) 2 


(A2) 


By evaluating the derivatives one obtains an explicit form 
for the Pauli-Jordan function in terms of the delta func¬ 
tions and their derivatives: 











(A3) 


A ij (x,x') = — 


3(Ax) * l (Ax )J 


- 5 lJ 


<5'(|Ax| — At) — <5'(|Ax| + At) 


|Ax| 2 

<5(| Ax | — At) — <5(|Ax| + At) 


|Ax| 2 


|Ax| : 


(Ax)’(Ax) J 

|Axp 


- 6 iJ 


<S"(|Ax| - At) - <S"(|Ax| + At) 


|Ax| 


where Ax = x — x' and A t = t — t'. The above deriva¬ 
tives of the delta function are to be understood as dis¬ 
tributional derivatives. An important remark is the fol¬ 
lowing. The expression we have obtained is ambiguous 
when At —> 0 or |Ax| —> 0 since they introduce prod¬ 
ucts of delta functions with functions which diverge for 
|Ax| = 0. In this case, one possible route is to consider 
the second line of expression (IA2|) with a finite e. 

Another important commutator of the electric field is 
the Hadamard’s elementary function. It is given by 


D ij (x,x') = (0\{E\x),E j (x')}\0) 

= 0 (1) 11 (t-*'.*-*'), (A4) 


dxt dx( 


where 


D^ l \t — t' ,X — x') = — —r-f 77-77-i-foT 

4tt 2 \ [(t — t' — *e) 2 — |x — x'| 2 ] 

1 

+ [(t — t' + ie ) 2 — |x — x'| 2 ] 

After calculating the derivatives one gets 

2(Ax) i (Axp - S ij [|Ax| 2 + (At - ie) 2 ] 
[(t — t' — ie ) 2 — |x — x'| 2 ] 3 
2(Ax) ! (Ax)3 - 6 ij [|Ax| 2 + (At + ie) 2 ] 

[(t — t' + ie ) 2 — |x — x' | 2 ] 3 

From the results above one may evaluate the Hadamard’s 
elementary functions found in equation m- On the 
other hand, from equation (IA3I) one can compute the 
Pauli-Jordan functions appearing in equation (fTdl) . 



D y (i, x') = — - 
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